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Abstract 


Nuclear interactions of high-energy alpha particles with target nu- 
clei important for cosmic ray studies are discussed. Models for elas- 
tic , quasi- elastic, and breakup reactions are presented and compared 
with experimental data. Energy- dependent interaction cross sections 
and secondary spectra are presented based on theoretical models and 
the limited experimental data base. 

1. Introduction 

In this paper we consider theoretical models of alpha-nucleus interactions for the purpose 
of developing a data base for describing the transport of cosmic ray alpha particles ( 4 He) and 
secondary light ions through bulk materials. Alpha (a) particles represent about 10—12 percent of 
the primary galactic cosmic ray flux, second only to the hydrogen component, which represents 
86-89 percent of the abundance. Interaction cross sections for 4 He-*H collisions have long 
been of interest (refs. 1-3) in the understanding of the chemical composition of the primary 
cosmic radiation because their interaction is the chief mechanism for 2 H and 3 He productions 
which are absent from the primary sources. More recently, the possibility of long-duration 
manned missions beyond the Earth’s orbit has resulted in an increased effort to understand the 
possible late and genetic effects of galactic cosmic radiation (ref. 4) and to accurately predict 
the exposures to be encountered. The large primary alpha component and also the important 
secondary alpha component from the fragmentation of heavier elements warrant a study of their 
physical interactions in spacecraft materials and tissue. 

The internal structure of 4 He, which is unique in several aspects, suggests an individual 
treatment of the physical interactions. First, this nucleus is the most compact with charge 
(matter) radius of 1.67 fm (1.33 fm) in comparison with the lighter nuclei 2 H and 3 He with 
charge (matter) radius of 2.1 fm (1.71 fm) and 1.88 fm (1.45 fm), respectively. The compactness 
of 4 He is an indication of the large binding energy at about 7 MeV/particle and results in a first 
excited state far removed from the ground state at about 20 MeV. Similar to the other light 
ions with mass number A < 5, no bound excited states occur for 4 He. We note that the first 
excitation level of all other nuclei occurs at only a few MeV. Nuclear correlations of dynamical 
content are expected to be enhanced because of this compactness which is, in part, due to the 
absence of Pauli exclusion effects and to the closed shell structure. The A = 2,3 nuclei are 
now amenable to exact solution for their internal wave functions using realistic nucleon-nucleon 
interactions (refs. 5-8). This is not true for 4 He, although variational Monte Carlo approaches 
have been used (refs. 9—12). Also, nuclear matter or independent particle models used for the 
study of heavy ions are expected to fail for 4 He where the concept of a mean field does not 
hold. Finally, the description of nuclear fragmentation for light nuclei will be in terms of a 
one-step process of direct reaction, rather than the abrasion-ablation picture used for describing 
heavy-ion fragmentation (ref. 13). 

The interaction cross sections for alpha projectiles consist of the elastic channel, compound 
nucleus reactions, stripping and pickup channels, and the fragmentation reactions. For cosmic 
ray studies these cross sections are needed over the energy range from a few MeV to above 
100A GeV. However, because of the dominance of atomic/molecular interactions at low energies 
(including the resulting short range of the ions) and the small primary flux above about 10A GeV, 
the energy range from above a few 10’s of A MeV to about 10A GeV is most important. The 
nuclear absorption cross section imposes an important constraint on the sum of the inelastic 
channel cross sections. A second-order solution (refs. 14 and 15) to the coupled channel equations 
(refs. 16 and 17) of the optical model in the eikonal approximation has been developed. The 
absorption cross sections and elastic scattering distributions for alpha-nucleus collisions are 



calculated by using a realistic model for the 4 He ground-state one- and two-body densities. 
The total scattering cross section is also evaluated. 

An important effect that results from the large separation energy for the first excited state 
in 4 He is a significant cross section for inelastic reactions with composite nuclei where the alpha 
particle remains intact; i.e., the target nucleus fragments while the projectile alpha particle does 
not fragment. This quasi-elastic type of reaction has been described by the optical model (refs. 18 
and 19) and is applied here to estimate the fractional contribution of quasi-elastic processes to 
the absorption cross section for common shielding materials. 

The fragmentation of light nuclei is often described by a single pole diagram (refs. 20-22) 
in which a single nucleon or cluster reacts on the target with the remaining piece of the ion 
acting only as a spectator. This type of impulse approximation fails when an extended region of 
the kinematical phase space is considered. We have developed an effective three-body multiple 
scattering approach (refs. 23 and 24) to describe the two-body dissociation of light nuclei on 
composite targets. The effects of elastic and inelastic fragmentation and final-state interactions 
are treated in this model. For 4 He this model describes the 3 H -p, 3 H e-n, and _d-d final states. 
Interactions of knocked-out clusters with the target are described by the optical model of multiple 
scattering, and thus the energy dependence of these 4 He fragmentation channels is described. 
We apply this model to calculate integrated cross sections and compare the calculations with 
available experimental data. Deuteron production cross sections from 3 He and 3 H projectiles 
and the breakup of 2 H may also be evaluated in this model. 

The two-body dissociation states and the quasi-elastic scattering account for roughly half of 
the absorption cross section when meson production is included. At low energies the stripping 
reactions are important, but their importance decreases at high energies because the amplitude 
involves the dissociation vertex in an energy-dependent manner. Extensive measurements have 
been made for pickup on *H for 4 He projectiles. We use the Serber model (refs. 25 and 26) to 
obtain a mass number dependent scaling of these experiments. Estimates of compound nucleus 
cross sections can be provided from the EVAP-4 code of reference 27 for light-ion projectiles 
below 100 A MeV. The remainder of the absorption cross section is shared by the a — ► dpn and 
a — ► npnp channels. Theoretical models of these reactions are hampered by the complexity of 
the many body final states and the fact that the vertex functions for these breakup modes have 
not been evaluated in any realistic model. We note that semiphenomenological methods have 
been used with some success for the a — > dpn breakup (ref. 28). 

The outline of this paper is as follows: First, the optical model is used to evaluate total, 
absorption, and elastic scattering cross sections that are compared with experimental data. 
The calculation of the quasi-elastic cross section for inelastically scattered alpha particles is 
then described, and predictions for common shielding materials are made. Second, the model 
of two-body dissociation of light ions is used to predict energy-dependent fragmentation cross 
sections. The Serber model for pickup and stripping is used to estimate these cross sections 
for alpha projectiles. Third, and finally, a survey of the available experimental data is made 
and combined with theoretical predictions to give energy-dependent parameterizations of cross 
sections for secondary 4 He, 3 He, 3 H, 2 H, and in alpha-nucleus collisions. A parameterization 
for the energy spectrum of secondaries is also presented. 

2. The Elastic Channel and Nuclear Absorption 

The evaluation of the nuclear absorption cross section proceeds from the elastic scattering 
amplitude and the optical theorem. In the Eikonal coupled channels (ECC) model (refs. 14 
and 16), the matrix of scattering amplitudes for all possible projectile-target transitions is given 

by 

7(q) = £zj d 2 b e* qb {e i *( b ) - l} (1) 


2 


1 1 nr ii m iii mu 1 1 i i in > uni iiiniiiiiiiiiu 



where barred quantities represent matrices and bold quantities represent vectors. Here, b is 
the impact parameter vector, q is the momentum transfer vector, and k is the projectile-target 
relative wave number. In equation (1), Z is an ordering operator for the z-coordinate which is 
necessary only when noncommuting two-body interactions are considered. The phase elements 
of X are defined by matrix elements of arbitrary projectile-target states of the operator 

/oc 

dzt aj {r a -Tj + x) (2) 

QJ X 

where /i is the nucleus-nucleus reduced mass, a and j label the projectile and target constituents, 
respectively, r is the internal coordinate, x is the relative coordinate with x = (b,z), and t a j 
is the free two-body scattering amplitude in the overall center-of-mass frame. For a projectile 
transition from quantum state n to n! and target transition from u to z/, we write 


y 


r+oo 

„'( b ) = 5Z ^ / dz <nv\t, 
ocj _00 


ClJ 
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where Ap and Aj denote the mass numbers of projectile and target, respectively. Equation (3) 
is written in terms of transition densities p as 


*ni/,nV( b ) = 2/fe X J dz J dTa dr J Puu'{ T j) PmA r a) *aj( r a + ~ *) 
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or in terms of transition form factors as 
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where F and G are the projectile and target one-body form factors, respectively. 

The two-body amplitudes must be related to their values in the nucleon-nucleon (NN) center- 
of-mass (CM) frame where the physical amplitude is determined by experiments. Making 
this transformation and noting that the z-integration in equation (5) can be performed formally 
if commuting interactions are assumed reduces equation (5) to 

*rw/,nV( b ) = 2 n-;- mN J ^ Q ^ ** ^ nn /(-q) G vi/( < l) /NN (?) ( 6 ) 


where /nn is the two-body scattering amplitude in the NN CM frame. Equation (6) is convenient 
for calculations since it is essentially a one-dimensional integration if the form factors are known. 

The second-order approximation to the elastic (EL) amplitude is obtained by including all 
transitions between the ground and excited states and assuming that transitions between excited 
states are negligible. Furthermore, the density of all excited (EXC) states is approximated by 
an average excited-state density. The phase matrix is then of the bordered form 


X(b) = 


*EL 

*00,01 

*00,10 

*00,11 

*01,00 

*EXC 

0 

0 

*10,00 

0 

*EXC 

0 

*11,00 

0 

0 

*EXC 


\ 
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( 7 ) 
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where = Xqo 00’ The charset eristic equation of this bordered, matrix is 

(X E XC - A)^ 0-2 (X E l - A)(X EX c - A) - T 2 ] =0 

where No is the order of X, A is the eigenvalue, and T 2 is defined by 


T 2 (b)= £ ^00, ^m^OO 

nor 


( 8 ) 
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The eigenvalues are then given by 


Al,2 = x(X E L + XeXc)± 


t 2 


2 (^EL ~ X E XC) 


1/2 


+ Y 2 
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with all others taking the value Xexc- The form of the eigenvalues allows us to treat the 
scattering system as an effective two-channel problem with 


X = 


X E L 


T Xexc 

Then, from employing Sylvester’s theorem we find that 


( 11 ) 


/cc( q ) - 2 7 r / 


e -iq b < exp 


*(x E l + Xexc) 


COS 


f r 7 r 2\ , v sin ( X DIF + T2 ) ^ 

( X DIF + T J + IXdIF (x2[)if + T 2)1/ 2 


- 1 


\d 2 b 


( 12 ) 


where the subscript CC denotes coupled channels and the difference (DIF) is given as 

Xdif = x( x el - *exc) 

An expansion of equation (12) reveals, as expected, that X E XC appears only in third-order and 
higher order terms in /\j\ (q) ■ As discussed in reference 14, a reasonable approximation to Xexc 
is to assume the ground-state density for the excited states. If X E XC ' s set equal to X E l we fi n< l 


A 2 ) 

Jcc 


,(q) ~ -7T- [ exp(-iq- b)[exp(iX E L)cosT - 1] d 2 b 
27 r J 


(13) 


The c oherent approximation (ref. 17) is recovered in the limit of small T. 

By using closure to perform the summations in equation (9), T 2 is given as 

T 2 (b) = ApAj f J d2( l d V e_!qb e "* q ' b /nn(q) /nn(q') 

X [-A P A T F^( q) q') G ( 1 ) (-q) G (1) (-q') + (A P - 1 )(A T - 1) F< 2 W) G (2) (-q, -q') 
+ (Ap — 1) FW( q + q') G (2) (-q,-q') + (Ap-l) F (2) (q,q') G (1) (-q, -q')] (14) 
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where F^ V) and (G^ and G^) are the projectile (target) one- and two-body, ground-state 
form factors, respectively. 

Townsend (ref. 29) has considered Pauli correlation effects between projectile and target 
nucleons. Here, the first-order elastic phase is written as 

*EL( b ) = *DIR( b ) “ *Ex( b ) ( 15 ) 


The direct (DIR) term is written as 

XmR(b) = / d2q e ' qb g(1) ^) /NN(q) (16) 


and the exchange (EX) term is written as 

* EX<b) “ I d2q eiqb f( ‘T q) G<1>(q) 

X p FJ d 2 q - „iq' b / NN ( q + q ') C(q') 


We use the parameterization of /nn 35 


(17) 


/nn(q) = — ^ ~ %N exp (-^9 2 ) 


(18) 


where is the relative wave number in the two-body system, a is the two-body scattering 
cross section, D is the slope parameter, and p is the ratio of the real part to the imaginary part 
of the forward two-body scattering amplitude. Values for the energy-dependent a, B y and p are 
found in reference 29. The correlation factor is found as 


C(q) = ~ e-’ 2 /" 2 (19) 

in reference 29 with d = 1.85 fm -1 . 

The total (TOT) cross section is found from the elastic amplitude by using the optical theorem 
as follows: 

47T 

°TOT = bn /(q = o) (20) 

Equations (13) and (20) show that 


roo r 1 

cttot = 47I y b db 1 1 - - exp[-Im(XEL + Y)] cos[Re(X E L + T)] 

~ \ exp[-Im(X EL - T)] cos[Re(X EL - T)]| (21) 


where Im and Re denote imaginary and real quantities, respectively. The total absorption (ABS) 
cross section is found by using 

°TOT = CT ABS + °EL ( 22 ) 

where <r EL is the total elastic cross section. Integrating equation (13) by using dfi « d 2 q/k 2 and 
equations (21) and (22) yields 

ijabs = 27t J b ^l — ^ exp(-2 Im X EE )[cosh(2 Im T) + cos(2 Re T)]| (23) 
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2*1. Model Form Factors 

The one-body form factor is written in terms of the charge (CH) form factor as 

F (1) (q) = F C n(q)/F p (q) (24) 

where F p is the proton form factor taken as e Tp <l ^ with r p — 0.86 fm. For 4 He, an excellent 
fit to the charge form is given by (ref. 30) 

F CH (q) = [1 - (0.316 9 ) 12 ]e-(°- 681 ^ 2 (25) 

The harmonic well model is often used for A < 20 where the charge form factor is 

Fch = (1 - sq 2 )e- a( ? (26) 

and values for parameters s and a are from reference 29. For nuclei where a Woods-Saxon 
density is appropriate (At > 20), 


'’ CH<r) = 1 + Jr-Ityc (27) 

An exact Fourier transform to obtain the charge form factor for a Woods-Saxon density may be 
found in a series solution (ref. 31) 


where 


47T 

^ch(<?) = — Po 4>{q) 


, x „ f — cos (Rq) 7 rc sin (Rq) Coth(7rcg) 

m = sinh(,r cq) + R 


_ 2c_ y'/_ j.m meg exp(-mR/c) 


} 


(28) 


(29) 


The series in equation (29) converges rapidly, and the first three or four terms are accurate for 
most applications. Values for the parameters c and R are taken from reference 29. 

The second-order calculations are difficult because the two-body form factors of the projectile 
and target must be known. We next consider the Jastrow method of correlated wave functions 
in order to model the form factors of 4 He. 


By using a coordinate system unconstrained by the nuclear center of mass, we introduce the 
model form factors Fm (which are related to the intrinsic form factors) given by 


F {1 \ q) = 


F {1 \ q,q') = 


Fff(q) 

(30) 

*cM(q) 

^M(q,q') 
^cM(q + q') 

(31) 
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with the harmonic oscillator CM correction assumed with 


Fcrvi(q) = exp 



( 32 ) 


where R is related to the oscillator parameter. 

The Jastrow method of correlated basis functions (ref. 32) introduces a correction factor to 
the wave function calculated in a single-particle potential model in order to take into account 
the effects of the short-range repulsive part of the nuclear potential on the wave function. The 
Jastrow-correlated wave function is written as (ref. 32) 


^( r l ■•■M) = ^(r-’-rA) ]J 

i>j — 1 


(33) 


where ^4 represents the Slater determinant for the ground state wave function and the 
correlation factor /(r^rj) is assumed to depend only on the relative separation of r* and Tj 
and obeys 


jo (|rj r j | — ► 0) j 

l 1 (l r * — r jl Laxge) J 


/( r»,rj) 

The two-particle density is given by 

p(r,r') = iV J |^(r,r', r 3 , • ■ • , r j 4 )| 2 rfr 3 • • • dr A 


(34) 


(35) 


where N is the normalization constant. 

The Jastrow correlation factor contains up to ^4-particle correlations. Since our considerations 
are for two-particle correlations, we consider a low-order approximation to this model (refs. 33 
and 34) and write the model two-body density as 


Pm(x,x') = N p„(x) p s (x.') | 5 (x,x ')| 2 
5 (x,x') = 


with 


(36) 


(37) 


where (3 will determine the correlation length. In equation (36), p s is a single-particle density 
assumed to be determined by the Slater determinant in equation (33). 

By using a single-particle wave function of the form 


’Mr) = y/N, s e~ r2 / 2R2 

the one- and two-body form factors of 4 He are found as (ref. 15) 


(38) 


= Eg i e-oVi-i 

i = 1 


(39) 


and 


Fjg^q') = e -9 2 / 4 ^ e -^ q' 


Z— 1 


(40) 
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with 


and 


where 


<3 

+ 

, jCN 

il 

(41) 

jx 

r-H 

i 

■tS> 

II 

(42) 

v i — a i (j^i / a i ^ 

(43) 

dj = fli/2aiVi 

(44) 

C\ = Ct / (aivi) 3 / 2 

(45) 

C 2 = —2Ct/{<12 v 2) 3 ^ 2 

(46) 

C 3 = Ct / i a 3 v 3) 2 / 2 

(47) 

C T = [(am) 3 / 2 - 2{a 2 v 2 r 3/2 + (w) -372 ] 1 

(48) 


2,2. Results for the Elastic Channel 

In figure 1 we show the charge form factor for 4 He with the experimental data from refer- 
ences 30 and 35. The solid line is obtained using the model form factor from equation (39) 
with l/R 2 - 0.65 fm~ 2 and /? — 2.0 fm -2 , and the dash-dot line comes from the parameteri- 
zation of equation (25). In figure 2 the two-particle density for 4 He is plotted against the relative 




Figure 1. Elastic charge form factor for 4 He. Figure 2. Two-particle density for 4 He as function of 

relative separation distance r. 


separation. The solid line is from equation (36), and for comparison a higher order Jastrow 
calculation from reference 34 is shown as the dash-dot line. Both models lead to about the 
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same healing distance which is directly related to the position of the minima in the one-body 
form factor (fig. 1) in our model. In the Jastrow model of reference 34, the CM constraint was 
neglected which may account for the differences in overall magnitude between the two models. 

In figure 3 the elastic cross section for 4 He- 4 He scattering at 635 A MeV is shown as a 
function of the invariant momentum transfer where t = —q 2 . The experimental data are 
from reference 36. The second-order model shows substantial improvement over the first-order 
model, especially at the second diffraction maxima where double scattering dominates and leads 


to excellent agreement with the experiment. In 
data from reference 37 at IT GeV. 



Figure 3. Invariant distribution versus momentum trans- 
fer (—t) for a+a elastic scattering at 635 A MeV. 


figure 4 a similar comparison is made with the 



Figure 4. Invariant distribution versus momentum trans- 
fer (— t) for a+a elastic scattering at IT GeV. 


Second-order calculations become more difficult for heavier targets if realistic form factors 
are used since many integrations must be handled numerically in evaluating T(b). In figure 5 we 
show a comparison between first- and second-order calculations for elastic 4 He on scattering 
at IT GeV with the Jastrow model described above used for 16 0 in evaluating T(b) using the 
appropriate radii for 16 0. Differences between the bordered and coherent model solutions are 
not substantial, which may be due to our choice of form factors for 16 0 in T(b). 

In table 1 the calculations of total and absorption cross sections are compared with experiment 
for several laboratory (LAB) energies (ref. 38) for a particles reacting with 1 H, 4 He, and 12 C. 
The second-order calculations are seen to lead to improved agreement with experiment; however, 
differences between the two solutions are only a few percent. We conclude that the second-order 
solutions offer improved agreement over the coherent approximation. However, this agreement 
is already within a few 7 percent for absorption cross sections. Improvements in elastic spectra 
occur only beyond the forward diffraction peak. Important improvements will likely be found 
at low energies. 
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Figure 5. Elastic 4 He+ 16 0 scattering at 1 A GeV. 


Table 1. Results for Total and Absorption Cross Sections 


Tlab> 
A MeV 

Values of err ot^ rob, for — 

Values of <tabs> rob, for — 

First 

order 

Second 

order 

Experiment 

First 

order 

Second 

order 

Experiment 

a-p 

870 

123 

140 

143 ± 1.6 

94 

101 

120 ±6.2 

2100 

126 

142 

147 ±0.4 

96 

103 

111 ±5.7 

aa 

870 

359 

389 

390 ± 6.3 

244 

253 

262 ± 18.5 

2100 

368 

397 

408 ± 5.5 

249 

259 

276 ± 15 

q- 12 C 

870 

814 

829 

790 ±7 

520 

528 

542 ± 16 

2100 

826 

842 

835 ±5 

530 

536 

547 ±3 


3. Quasi-Elastic Scattering 

The large energy separation between the ^He ground state and first excited state may lead 
to a large cross section for quasi-elastic scattering. In the quasi-elastic process the alpha 
projectile loses energy and receives momentum transfer without suffering a change in mass and, 
concomitantly, the target nucleus fragments. A significant quasi-elastic (QE) cross section will 
be important since this process will contribute to the absorption; however, no secondary particles 
are produced from the projectile in the reaction. The quasi-elastic distributions in momentum 
and energy transfer are also used to describe the inclusive breakup of light ions (refs. 23 and 24). 
The QE cross section is evaluated next in the high-energy optical model. 


3.1 Scattering Formalism 

In treating inelastic scattering we assume that the off-diagonal terms in X, denoted by X 0 , are 
small compared with the diagonal one, X D , and then we expand / from equation (1) in powers 


of X 0 ' 


-yJ 


d?b e 


i<l'b e iX D {b) ^ 

771=1 


[i X 0 (b)]’ 


ml 


} 


(49) 
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We also will make the assumption that the diagonal terms are all represented by the ground- 
state elastic phase X. In the remainder of the paper we drop the subscript EL on the elastic 
phase because it will appear only in a distorted wave, and its meaning should be apparent to 
the reader. By using equation (3) we sum over target final states X (continuum) to find the 
inclusive (IN) angular distribution for the projectile when its mass remains unchanged as 

= (S? J dH dh ' eiq(b ” b,) exp,i[ * <b) - *>'>]> 

x 5Z S |n P .Y > 

I/0m=l Km -> 

x < wY0p| iX"*"(b')j |0p0y > (50) 

Equation (50) allows only for a study of the momentum transfer spectra of the projectile. In 
considering the energy loss of the projectile, energy conservation must be treated. By using 
continuum states for the target final state, energy conservation leads to 


d 2 o \ 

dQ dEp 


(27T)- 


r A T 

I d 2 b d 2 b f e iq '( b ~ b ) exp {i [X(b) -X+(b')]} ^ W m ( b, b 7 ,^) 

m— 1 


(51) 


where Ep is the energy of the projectile in the final state and uj is the energy loss of the projectile. 
We define 


i f r die 

w m (b,b', 0 j) = n [p^J HE, ~ Ei) < OpO r |p(b)] m |o p k J - 


> 


i TTl 

x < kjOp X (b ) OpOp 


> 


(52) 


(refs. 18 and 19) where kj is the wave number vector of a knocked-out target nucleon. The 
functions W m are next related to the response functions of the target in the cylindrical geometry 
of the eikonal approximation. 

The first collision term is written as 


tU(b,b» = J d 2 q d 2 q' e^ b e~^ h ' F(q) F(q') 

x /nn(q) /nn^ ) J ?27rj2 ~ ^k) ^ok(q) G+(q') ( 53 ) 

where F is the projectile ground-state form factor and is the target transition form factor. 

We change variables as 

« = ^(q + q') (54) 

0 = q - q' (55) 

x = s - s' (56) 

y = ^(s + s') (57) 
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and also 


such that 


R = b - b' 

S = i(b + b') 


(58) 

(59) 


Wj(R, S,cj) = J d2a d2 P e * a R e?/3 S A ( a + f ) A+ ( a “ f) ( 60 ) 

A(q) = F(q) / NN (q) (61) 


where we have defined 

and the target response function is 


Rl(a,0,u) = J j^-2 S(u - F k ) G ok ^ G£ 0 - 0 (62) 

By following Krimm et al. (ref. 39) we can formally treat the delta function in equation (62) by 
introducing a Fourier transform pair 


Then, 

dk 

(2tt) 2 

For a nonrelativistic nucleon we have 


Rl(a,(3,t) = J 


= J ^e iut Ri( a> /3,t) 


(63) 

J dw e~ iut R\ (a,/3,uj) 


(64) 

~ iEk ‘ G„ k (a + G+ | 

M) 

(65) 

k 2 

" 2m jv + CBl 


(66) 


where is the binding energy. By assuming plane waves for the target final state in G ok , 
equation (65) then becomes 


RfaM = J J^odx dy e- ieB i l e ikx 


( 2^) 5 


x $ 


y +x U+ 


M 


(67) 


where $ is the single-particle wave function of the target ground state. Using equations (67) 
and (63) gives 


' ^ fdx dy e ia x e i0 y - e Bl )x 2 

Hl( a tA w ) = < x $ (y + f ) 4> + (y — f ) (w > ejjj) 

0 ( w < e Bi) J 


> ( 68 ) 


The higher order collision term is more complicated because of the enumeration of projectile- 
target intermediate states that can occur. A first approximation is to keep only Ip-lh excitations 


| 


■ 

I 
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of the target (one for each inelastic scattering) and assume that the projectile remains in the 
ground state (coherent approximation). 

Using similar coordinate changes as described above, the rath-order collision term is found 
in the coherent approximation to be 


a 2m Am r m r 

- (m i) 4^ / n ** 

X R e l,3 r s Aj (aj + A+ (etj - 
x i, . . . , (3 \ , . . . , /3 m , a;) (69) 


where 


ra m /* m . 

Rm ( a 1 j • ■ • > a m ? /^1 , . • • > /3m > = J TI ^ X J ^ % e J 


■ x j 


x * (y> + y) * + (yj - y) 

om— 1 / ^ \m— 1 

2 


X 


2 myv ( u > - efl TO ) E 

J=1 


-l)/2 


X 1 


I 2 mN - e Bm ) ^ 


i=i 


(70) 


where i? m = 0 for u) < e Bm . We next consider a simplified representation of the ra > 1 terms. 
By assuming that the target wave functions are forward peaked, we approximate 


Jm— 1 y m y ^ ^jy 

/ 1 \ 771-1 




(m- 1)! 2 m_1 . -*■ 

7 = 1 


-1 II J ° o(m-l)/2 + ° (^*®i) ( 71 ) 


2( m -!)/ 2 


where 
such that 


= yj 2tojv - f D m ) 


/ \m— 1 771 

Rm{°^ h • - • j %) /3i, . . . , f3 m ,uj) ~ IT 


and 


(ra — 1)! 


/ _ \m-l , 

w m (R, s,w) = y £ f™/ 2 ^ wi 

(ra — l)!(ra!)^ 


3 = 1 


R, S, 


^71 


a j'@j' 2 (m— 1)/2 


2( m -!)/2 


]}’ 


(72) 


(73) 


(74) 
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A numerical test of the forward-peaked wave function approximation is discussed below. The 
energy loss spectrum in a coherent projectile model is given by 


d‘ 


dn dE 


= -^2 f d 2 R d 2 S e* q R exp{ i [ X (^) _X+ ( 
lE'p) in (2tt)2 J 


S-R 


)]} 


At / \ m— 1 

x ST ^~ eg m) J w 

4 (m - l)!(m!) 2 ^ 1 

m = 1 


{' 


R, S, 


2 ( m - 1 )/2 


]}' 


(75) 


After angular integration, the energy loss spectrum is found as 

$) m - i* s ,m x(si i few 4 «] r <7e) 

The coherent approximation assumes that the projectile remains in the ground state through- 
out the scattering. The leading-order correction to the coherent terms occurs in W 2 and corre- 
sponds to the replacement (ref. 40) 

4 F («i + y) F ( Q i - t) f (“2 + y) f (“2 - y) 

- Ap | F( 2ai) + (A P - 1) F (a : + y) ^ («1 - y)] 

x F(2q 2 ) + (Ap - 1) F (a 2 + y^ ^ ^«2 - y^J j ( 77 ) 

which, physically, represents the projectile dissociating in the intermediate state. Further 
modifications are necessary when correlation effects not included here are treated. 

The distribution in momentum transfer to the projectile nucleus may be obtained from 
equation (75) after integration over the energy loss. An alternate expression is obtained from 
equation (49) or (1) through use of closure without regard to energy thresholds for ejecting 
particles into the continuum given by references 40 and 41. Thus, 

£)„ - (s)7 * R Ss e ‘ qR exp [ x ( s + 1) - x+ ( s - !)] } [ e!><RS) - *1 (78) 

with 


fi(R, S) = t ^|^2 / d ^ d V e * qb e_iqb '/NN(q) 4 N (q')^(q) F(q') G{ q - q') ( 79 ) 

(27 rltNMf J 
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The distributions in momentum transfer obtained from equations (75) and (78) are equivalent 
if the response function obeys the sum rule 

J du J?i(q,q',w) = G(q-q') (80) 

with similar relationships for higher order terms. 

An accurate approximation for obtaining numerical results is to expand the elastic coupling 
phase in equation (75) or (78) as 



= —2 Im X(S) — zR * Vg Re X(S) + . . . 


(81) 


with 

zR • Vg X(S) = — Rcos(<f)ji — 4>s ) R e X^S) 
where <fij{ and (f>g are azimuthal angles and 


X'(S) = [°°q 2 dq MqS ) / NN (g) F(q) G{q) 

*NN Jo 

We then have, for example, 


<fi<7 \ 

di 1 dEp J 


ju 2 roc 

-—2 / RdR 
(2 *T JO 



S dS e ~ 2 Im X (5) J 0 ( q R)J 0 {R [R e x'(S)]} 


^ (m — l)!(m!) 2 

771— 1 / v / 



R) Ri 


(m — l)/2 


]}’ 


(82) 


(83) 


(84) 


The second term in equation (81) physically allows for momentum transfer in elastic scattering 
and usually makes only a small contribution (ref. 41). 


3.2. Shell Model Response Functions 

For light nuclei (A < 16), we use shell model harmonic oscillator wave functions. Thus, for 
s-shell nucleons, 

/ \ 3/4 

** (r > = (^) e_r2/24 (85) 

and for p-shell nucleons in a spherical basis with components m, 


4Wr) = 

where is the target radii and 


.) 


3/4 


ttRj, J V Rj, 


rm e-^ 


r 77l ~ 


- / ^( x±iy ) = 


*"} 

= 0)J 


t z (m = 

The s-shell and p-shell probabilities are given, respectively, by 

C 4 


( 86 ) 


(87) 


( 88 ) 
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and 


(89) 


Cp — 


A- 4 

A 


From equation (68) the response function is then 


Rl(a,/3,^) = m N Rq- e _jR T^ /4 e ^ e R T^t 

x | 2 C s - + \ C P Rt (« 2 + *l)] ^ ( 2 4“?l) 

- |cp/4«6 A (24^1)} 


(90) 


where I Q and Ii are modified Bessel functions. Higher order response functions are then 
approximated by using equations (73) and (90). The collision term W\ can now be found 
in analytic fashion, and higher order terms are approximated by using equation (74). 


3.3 Results for Quasi-Elastic Scattering 

Theoretical calculations for quasi-elastic a 4 He scattering at 1.4 GeV are compared with 
experiment (ref. 42) in figure 6 for several scattering angles. The calculations shown are 
made with the approximation of equation (75). For consistency we show only the coherent 
contributions since we have not formulated the corrections for incoherent projectile motion 
beyond W 2 ■ The multiple scattering structure is apparent with single inelastic collisions dom- 
inating at a small momentum transfer and the higher order contributions increasing in im- 
portance with q. In figure 6 the second collision term is seen to peak at a smaller energy 


I 


First collision term 

Second collision term 

Sum of first and second terms 

— Sum through third collision term 

° Experimental (ref. 42) 





B 

I 


i 


Figure 6. Momentum spectra of a particles in a— 4 He collisions at 14 GeV for scattering angles of 2.112° 
(q = 1.31 fm -1 ), 3.094° (9 = 1.92 fm' 1 ), 3.63° (q = 2.25 fnT 1 ), and 4.552° (q = 2.82 fm” 1 ). 
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loss than the first term because of the sharing of the total momentum transfer between two 
scatterings. This effect keeps the position of the theoretical peak in good agreement with the 
experiment. The strengths of the distributions compare fairly well with experiment except at 
the largest momentum transfers. This discrepancy grows when calculations are compared with 
the larger angle data of reference 40 (not shown), and it is attributed, at least partially, to our 
use of a Gaussian wave function for the 4 He ground state. 

In figure 7 we illustrate the accuracy of the approximation of equation (73) for the term W 2 
for a-a scattering at 0 = 3.63°. Incoherent effects are shown to significantly reduce the second 
collision term. Correlations among projectile nucleons will thus play a role in understanding the 
quasi-elastic peak. 


Coherent model with 
forward -peaked 
approximation (eq. (73)) 



Figure 7. Comparison of calculations of W2 in various approximations for a- 4 He scattering at 6 — 3.63°. Solid curve 
is exact incoherent result, dash-dot curve is incoherent result with forward- peaked wave function approximation, 
dashed-line curve is exact coherent result, and dotted curve is coherent result with forward- peaked wave function 
approximation. 


In figure 8 we show results for quasi-elastic alpha scattering on 16 0 at 1 A GeV for scattering 
angles of 1° and 4°. The higher order terms are seen to grow in importance with increasing 
energy loss and momentum transfers. 

In figures 9-11 we compare our calculations with the experiments of Ableev et al. (refs. 43 
and 44) for the inclusive scattering cross section versus the invariant momentum transfer on 12 C, 
27 A1, and ^ 4 Cu targets at 3.6A GeV. The measurements correspond with the sum of the elastic 
and quasi-elastic cross sections. In the figures the calculations are denoted as elastic scattering 
using the first-order optical model (dotted line), the quasi-elastic part using equation (78) 
(dashed line), and total scattering which is the sum of the elastic and quasi-elastic cross sections 
(solid line). Agreement with experiment is seen to be quite good. 

In figure 12 we show predictions for the total quasi-elastic cross section for a+ 12 C scattering 
as a function of laboratory energy. The energy dependence of the cross section follows roughly 
that of the nuclear absorption cross section and is seen to make up almost 10 percent of the 
absorption. In figure 13 we show predictions for a+ 16 0 scattering where similar conclusions 
apply 
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Collision terms 


First 

First and second 

First through third 

First through fourth 




(a) 9 = 1°. (b) 0 = 4°. 


Figure 8. Calculations of quasi-elastic a- le O scattering at 1 A GeV for 6=1° and 4°. 




Figure 9. Inclusive 4 He+ 12 C scattering distributions at Figure 10. Inclusive 4 He+ 27 Al scattering distributions 
3 6/1 GeV versus invariant momentum transfer. at 3.64 GeV versus invariant momentum transfer. 
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Figure 11. Inclusive 4 He+Cu scattering distributions at 
3. 6 A GeV versus invariant momentum transfer. 



Figure 12. Prediction of quasi-elastic c*+ 12 C total cross 
section versus laboratory energy. 



Figure 13. Prediction of quasi-elastic o+ 16 0 total cross section versus laboratory energy. 


4. Light-Ion Fragmentation on Nuclear Targets 

The fragmentation of 4 He is somewhat simpler than that of heavier nuclei in that there are 
only a small number of final states that can occur. These reactions are 


r3 He + n + I (91a) 

3 H + p + X (91b) 

a + T —> d + d+X (91c) 

d + n + p + X (91d) 


(n + n + p + p + I (91e) 


where X is the final target state. Each of the reactions in equations (91) can occur with or 
without meson production if sufficient energy is available. The reactions in equations (91) are 
not exhaustive of the absorption processes; most notable are the compound nuclear and pickup 
channels that are important at low energies. A model for the two-body dissociation of light ions 
has been developed (refs. 23, 24, and 45) that describes the first three reactions in equations (91), 
which we now discuss. 
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4.1. Two-Body Dissociation Theory 

For an inclusive reaction involving the two-body dissociation of the projectile, we write 


P + T-*a + b + X (92) 

where a and b are assumed to be clusters present initially in the projectile and X is the final 
unobserved target state. We consider the case where a is the observed projectile fragment in the 
measurement and note that the unobserved target states must be summed over in evaluating 
the cross section. A summation over possible states of the particle b should also be considered. 

By using relativistic kinematics, the transition amplitude (' T ) for equation (92) is related to 
the Lorentz invariant momentum distribution for producing the fragment a by 


Ea 


dc 7 

dp a 


= E ( 


( 27 tY 
1 0 


/ A? m 

dpb eii dp j S( pj - Pt) 6(Ef - Ei ) \TjiY 


(93) 


where / and i label the final and initial states, respectively, (3 is the relative projectile-target 
velocity, and the summation in equation (93) is over the possible configurations of target particles 
in the final state. A useful approximation to equation (93) is to consider the final target state 
as an effective particle X and apply energy conservation in an approximate manner. Here we 
consider an effective three-body problem 


where the phase space factor is given by 

„ EgEjjExpl 

~ p b (E b + E x ) + PaE b cos (6 a + 0 b ) 


(94) 


(95) 


In equation (94) we are ignoring the mass spectrum of the final target states in imposing energy 
conservation, and we will assume that M x ~ Mp. 

The transition amplitude can be written as a three-body problem of a — T, b — T, and a — b 
interactions when rearrangement channels are neglected and with the understanding that all 
target final and intermediate states must be summed. Using the Faddeev method allows us to 
consider the multiple scattering series generated by the coupled set of integral equations 

T = T a + T b + f T (96) 


T a = % T + f bT G 0 (? b + T t ) 
f 6 = f aT + f aT G 0 (T a + T t ) 
f T = f ab + f ab G 0 (f a +T b ) 


where f aT , f bT , and T ab are the “two-body” amplitudes that are the transition operators for 
aT, bT, or ab scattering, respectively, in the projectile target space and where the Green’s 
function in the impulse approximation is 


G 0 



k 2 

K b 

2 m b 



1 


(97) 
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We consider the leading-order corrections to the pole approximation by truncating equation (96) 
as 

T= (l + T ab G c 3 ) (f aT + f bT + f aT G 0 %T + %t G oTarj (98) 

and replacing T a x an d ^bT by their on-shell values. However, equation (98) allows for all orders 
of multiple scattering by assuming the dominance of the ab cluster in the projectile and the fact 
that ab scatters only after interaction with the target. 

The first-order terms of are shown in figure 14. In figure 14(a) the fragment a is the 
spectator with the unobserved fragment b interacting with the target. In figure 14(b) the roles 
of a and b are reversed with b being the spectator. These terms are written as 

Tjp = <f>( u a ) T bX (^, Q) + <£( n b ) T aX (y/stir, Q) (99) 

where (f) is the overlap function representing the virtual projectile decay, P —> a + b, y/s is the 
invariant energy for the quasi-scattering of cluster on the target, Q is the total momentum 
transfer in the reaction Q = — pjy, and 


U a = p a P p (100) 

rri p 

u b = ~Pb + — Pp = Q - Pa (101) 

nrip 

From equations (100) and (101) we expect the first term in equation (99) to dominate at small 
p a . In equation (99) the amplitudes Tj T are half-off shell, where j is either a or 6, but are 
assumed on-shell in a high-energy approximation. 



(a) Spectator term. (b) Participant term. 

Figure 14. Terms for projectile fragmentation. 


In figure 15 corrections for final-state interactions (FSI) between the projectile fragments are 
shown. The FSI diagram (fig. 15(a)) leads to the integral 





4> (k + ^Q) T ab { k, p ab ) 
p lb - k2 + ie 


( 102 ) 


where p,,/, is the relative momentum. We follow references 46-48 and use an off-shell separable 
T-amplitude for the FSI where the overlap function is used to replace the Yamaguchi poten- 
tial form factors, thus ensuring orthogonality between the bound and scattering states of a and b. 
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(a) Spectator term. 


(b) Participant term. 


Figure 15 . Terms for final-state interaction. 

The higher partial waves are estimated using the Glauber amplitude for ab elastic scattering. 
The terms in figures 14 and 15 are combined as 


TfP — 4>(u a) T bx {y/svr, Q) + <t>( u b) T a x(V s aT> Q) 


(103) 


where the distorted overlap function (p is defined as 

D (p ab , 2 -^q) 

<t>(u a ) = 4>( u a) - 4>(Pab)' ~ 


D(Pab> O) 

- J sin 9 cos 9 dd <p (p ab k + q) f a b (%Pab sin ^ ) 


and 


Dip, q) — ~ J 


dk 


<t> (k + S[) (q 2 + k 2 ) </>(k) 


(104) 


(105) 


p 2 - k 2 + ip 

where a is related to the a - b separation energy e s by a 2 = -2p ab e s . The energy-dependent 
parameters of the two-body amplitudes in f ab are taken from reference 29 with the two-body 
cross sections allowed to smoothly fall to 0 below 20 MeV such that the s-wave part dominates 
at low energies. 

In equation (103) the decay amplitude <f> is written in the CM frame; however this amplitude 
is usually parameterized in the projectile rest frame. We transform to the projectile rest frame 

^ I/O 

, , \ — ^aiEp ~ 

^ Uaj " [E' a (m p -E' a )m p \ 

where primed variables represent projectile frame quantities. Using the parameterization 


<f(Pa) 


(106) 


l—\ l 


(107) 


with a\ — a and V~N being a normalization constant allows the dispersion integral in equa- 
tion (105) to be evaluated in analytic form. Values for the vertex function parameters are listed 
in table 2. The overlap probability |Z| 2 is also listed in table 2. 

In figure 16 we show the rescattering-type correction where both a and b interact with the 
target. In figure 17 the corrections for FSI to the rescattering are shown. These terms are almost 
always neglected for composite- composite breakup reactions. For figure 16 we write 


rj? = 2/i aT 


f , TfjpiQ - q 2 ) 4 > (pa — ^Pp Q) 7 1 a 7 ’(q 2 ) 

J d< l2 ^2 


-Q 2 - ( 2 P aX ■ 02 ) + tn 


(108) 
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Table 2. Overlap Function Parameters 


Projectile 

E s , MeV 

a\ y fm 1 

8 

r 

1 

i 

03 , fm 1 
(a) 

|Z | 2 

d — ► 7i - p 

2.22 

0.232 

1.434 

0 

1.00 


7.10 

.448 

.92 

0 

.85 

3 He —> d — p 

7.10 

.420 

.92 

0 

.85 

4 He — » 3 H — p 

19.82 

.846 

1.15 

1.65 

.60 

4 He -» 4 He - n 

20.58 

.863 

1.65 

1.65 

.75 

4 He d — d 

23.85 

1.070 

1.60 

2.50 

.50 


“For a 3 = 0: a \ = 1; a<i — —1. For a 3 ^ 0: a\ = 1; — 



03 = —a 1 — a 2 - 



Figure 16. Rescattering corrections. 



(a) Spectator term. (b) Participant term. 

Figure 17. Rescattering corrections with final-state interaction. 


with a similar contribution occurring when the roles of a and b are reversed and where p a x is 
the relative a and X momentum in the final state. 

For figure 17 we find that 


r (2)/ _ 



TabQc-abi Pafr) 

p lb- k lb + ie 


X 


/ 


dq 2 


TyAQ - q 2 ) 4> (Kb - q 2 - go) ^ar(q 2 ) 

2^1 - ' (fir - Si )] + ka6)2 + i7] 


(109) 
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Equation (109) is considered so that the rescattering corrections can be treated in a consistent 
manner with the first-order terms because the fragment FSI will occur at a relatively low energy 
compared with the aT or bT motion. The contribution from equation (109) should be most 
important for small p (J ^, and we assume that p\ P a b or k ab ' n order to reduce the integral to 
a manageable form. We can then combine equations (108) and (109) in a distorted-wave form 
as 

?}?**-— ^ /xdx%(Q-x)?(x,p a ,Q)T ar (x) (110) 

Jl PaX J 

where 


<f> (x, Pa, Q) = 4 > ( Pa - — P p - X ) 
V m p / 


2Pab PPTPaX 

PaTPX 


X 


/ 



Pab) 

Plb ~ k2 + iV 


( 111 ) 


where x = —2 p a x cos 0 q2 fa, and the second term in equation (111) is evaluated similar to that 
in equation (104). In equation (110) the target final and intermediate states must be considered 
for inelastic fragmentation. 


4.2. Elastic and Inelastic Fragmentation 

A convenient way to handle the target state summations is to separate the momentum 
distribution in equation (94) into elastic and inelastic terms corresponding to the final target 
state. In elastic fragmentation the target remains in the ground state, and in inelastic 
fragmentation it remains in the target fragments. This separation is given as 



( 112 ) 


Using 

fa = ^ T a ( 113 ) 

where ( XJ — E x Ej /(E{ + Ej) and defining the internal momentum distribution 4> of the fragments 

(in units of (MeV /c) “ where c denotes the speed of light) by <fi = (2tt )^/ 2 T allows us to write 

for elastic breakup 


E a 




1 

P 


/ 


dSlp 1C |M EIj | 2 


(114) 


where 


Mel = — $( u a) fvr( Q) — — $( u f>) /or(Q) 

*bX e aX 

+ o — [ x dx /«r(Q - x)$(x,p a ,Q) / oT (x) 

2 taX^TPaX J 

+ o — — — [ v d v /ot(Q - y)^(y>P6-Q) /tr(y) ( 115 ) 

2 t a T%XPbX J 


where f aT is the elastic amplitude that is evaluated in the coherent approximation to the optical 
model. 


E 

K 

i 

i 
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For inelastic breakup we have 


E a 



( 116 ) 


where the first-order terms are given as 


Ml 2 - M 2 *L(Q) + l^(«6)| 2 ^-(Q) 


El M ml 2 = ^l^ u « 

X^O C bX 


d£l in 


aX 


dSl Ijsj 


+ 


^ {*(.!„) S+(u s ) +<!.+ («,,) $(«,,)} ^(Q) (U7) 


e aX e bX 


and with the definitions 


da a 

dfijpvj 


^|<Xa|/ ar (Q)|Ta >| 2 


(118) 


and 


da ah 
dill n 


E < aT\f+ T maX X Xb\J bT (Q)\bT > 

x/o 


(119) 


Equation (118) is just the inclusive distribution for the reaction a + T — ► a + X which 
corresponds to equation (78) when the optical model is used. The cross section of equation (119) 
is an interference effect that occurs when there are two virtual projectile clusters available 
to fragment the target. At high energies we use an on-shell approximation to evaluate 
equation (119) which is an extension of the optical model result in equation (78). Thus, 


_ PaXPbX f d 2 b d 2 b / iq (b-b') exp { i [^ar(b)-X+ T (b')] } Lf2 a6T (b,b') _ J (120) 

dS \ N (2t r) 2 J F L J v ; 


where 


n abT = ^^2 J d ^ d V eiq b e_iq ' b ' /NN(q)/^ N (q')^a(q)F 6 (q / ) G(q - q') (121) 


The higher order inelastic terms are numerous and include terms where both a and b excite 
the target, only a single cluster excites the target with the second scattering elastically, and 
the interference terms between the second-order and first-order terms that lead to identical final 
target states. These terms must include both orderings for a scattering prior to 6, and vice versa. 


25 



As an example of such processes, double scattering in which only one projectile cluster excites 
the target nucleus is written as 


2 

lAfffl 2 = (- — — J [ x dx x dx' $(x,p a ,Q) $ + (x',p a , Q) 

1N \^aX%TVaXJ J 


*'> *'> + ^ ~ * ' * ' It ' ) 


2 

+ ( - — — J [ y dy y'dy' $ (y, p&, Q)$ + (y', Pb, Q) 

\^ e aT e bXPbX / J 


’O"') + ~ yQ - y,> 


where we have defined 


dcr 

d^l el 


(x, x') = /(x) / + (x') 


(122) 


(123) 


and 


dcri 

dQ i^sT 


(x, x') = / dH d2b ' eiX b e_<X ' **' exp b [ X ^( b ) “ X J + r( b ')] } (e°^ (b ’ b,) - l) 

( *' J (124) 


4.3. Impulse Terms for Inclusive Deuteron Production 

In an inclusive measurement the reactions in equations (91c) and (91d) are not distinguished. 
The reaction in equation (91d) is more difficult because of its four-body final state with at least 
three relative motions needed to be considered. Also, the vertex function for *11 e * npd has not 
been evaluated. We estimate the energy spectra for equation (91b) by considering the impulse 
terms for this reaction. Higher order scattering terms contribute largely as a normalization 
correction, and effecting the tails of these spectra will be considered elsewhere. 

In figure 18 we show the impulse terms for inclusive deuteron production. These terms are 
written as 


Tfi - <!>d(Pd) T d X ( Q) + <£d(Pd _ Q) T dr( Q) + <t>dn{Pd,Pdx) T px( Q) 

+ <t>dp(Pd’Pdx) T nX(Q) + 4>np{Pd,Pnp) T dX (Q) ( 125 ) 

where fa is the a — > dd vertex and fa is the a -* dnp vertex. We approximate the 4>%j vertex 
function by assuming that the weak binding of the deuteron is such that a correlated neutron- 
proton pair in 4 He closely resembles a deuteron cluster in the 4 He ground state. We then assume 
that 

J $d P (.Pd, P dp) dpdp w fa(P d) ( 126 ) 

The calculation of the momentum distribution for inclusive deuteron production then closely 
resembles the evaluation of the impulse terms for the two-body dissociation. We note that the 
reactions in equations (91c) and (91d) are orthogonal in the projectile Hilbert space. 
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Figure 18. Impulse diagrams for inclusive deuteron production. 


4.4. Results for ^He Breakup 

We first compare our results with experimental data for 1 H targets where inelastic breakup 
does not occur. In figure 19 the angular distribution for 3 He production in 1.02 A GeV a particles 
on ! H collisions is compared with the experimental data of Bizard et al. (ref. 49). Contributions 
from various breakup terms are labeled on the figure with the total (solid line) agreeing quite 
well with experiment except at the largest angles. 

In table 3 predictions of total production cross sections for 3 He and 3 H are compared with 
experimental data from references 49—51. The data from Webber (ref. 50) are preliminary. 
We note that Webber reports that a p < ct.jjj; this cannot be correct because the mechanism 

for producing 3 H will always produce a proton, and several other mechanisms for producing 
protons exist from equations (91). An error in distinguishing Z = 1 fragments is probably at 
fault. Agreement between theory and experiment is satisfactory. Calculations were made with 
fixed-energy NN parameters. This ambiguity in on-shell amplitudes will be most important at 
lower energies (<500 MeV) and for light targets where slope parameters are more important, and 
thus it should be studied in more detail. The absence of pion production in the model prevents 
a realistic comparison much higher than 1 A GeV for total cross sections. In table 4 predictions 
are given for the production of 3 He from ^He and for deuterons from 3 He on several targets. 
Vertex function parameters for 3 He -> pd are estimated from Kok and Rinat (ref. 52). For 
composite targets the higher order terms in equation (98) were not included because of the large 
computational time required. The higher order terms are expected to increase in importance for 
heavier targets and lower energies. 

Calculations of the longitudinal momentum distribution for triton production at 1.9 A GeV 
are compared with experiment (ref. 22) in figure 20. The dash-dot line is the plane-wave- 
impulse approximation (PWIA) for proton knockout and clearly underestimates the data. 
This discrepancy would be only partially resolved by using a wave function with higher 
momentum components. The short dotted line and the dash-dot line are the impulse terms 
with FSI for proton and triton exchange, respectively. Note that the final-state interaction 
causes enhancement in the cross section at large momenta. This conclusion was also found 
in reference 53 which used a Gaussian wave function and the Glauber model (ref. 54) for 
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spectator 

participant 



Figure 19. Angular distribution for a-f- 1 !! — » 3 He at 1.02 A GeV. 


Table 3. Comparisons of Calculations With Experiments of Webber (ref. 50) 
for A = 3 Fragment Production From 4 He . 

[Calculations are given in parentheses] 

(a) a+ 12 C — + Ap (b) a+hl ^Ap 


T L ab, A MeV 

<73 He , mb 

377.1 

26.3 ±2.6 (19.5) 

519.9 

26.4 ±2.6 (20.8) 

1025 

24.1 ±1.9 (22.5) 


Tlab, A MeV 

03 Hl mb 

<T3 He . mb 

203.3 

93.1 ±9.3 (77.2) 

60.4 ±6.0 (79.3) 

377.1 

79 ±7.9 (59.9) 

66.9 ±6.7 (60.9) 

519.9 

(62.1) 

69.4 ±6.9 (59.8) 


Table 4. Calculations of Fragmentation Cross Sections 
for Light-Ion Breakup 


(a) a+A T -* 3 He ( b ) 3 He±A T — 2 H 



Values of (72 H , mb, for — 

Tlab s A MeV 

16q 

27 A1 

500 

158.0 

211.2 

1000 

128.5 

184.0 



Values of <X3 He > mb, for — 

Tlab. A MeV 

12 C 

16 0 

27 A1 

500 

59.8 

72.3 

96.8 

1000 

64.1 

74.7 

100.1 
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1.05 A GeV a particles. The off-shell s - wave part of the FSI dominates this effect. Including the 
Glauber amplitude for higher partial waves represents only about a 10-percent correction. The 
inflection in the dash-dot line near 0.2 GeV /c is caused by the interference between the s-wave 
and Glauber terms. The solid line and the dashed line in figure 20 are the sum of all terms with 
and without the effects of interference between the diagrams, respectively. The full calculation 
provides a good description of the data out to 0.4 GeV/c, but it underestimates at larger values. 
The effects of interference between diagrams are clearly important. 

In figure 21, calculations of the transverse momentum distribution for a beam of energy of 
2.09A GeV are compared with the experiment of Anderson et al. (refs. 55 and 56). Triton 
exchange is seen to contribute only for small values of pp. The proton knockout term with 
FSI provides a good representation of the data out to 0.6 GeV/c and falls below the data at 
higher values. Calculations of pion production in deuteron breakup (ref. 57) suggest that we 
should expect some contribution from this mechanism at large momentum. In figures 22 and 23, 
calculations of the transverse momentum distribution at 0.385 and 1.041 A GeV, respectively, 
are shown in comparison with data from references 56 and 57. The total calculation agrees well 
with the experiment. 


PWTA (proton exchange) 



p, GeV/c p f GeV/r 

Figure 20. Comparison of calculations of longitudi- Figure 21. Comparison of calculations of transverse 
nal momentum distribution for 3 H production in momentum distribution for 3 H production in a— 12 C 

a— 12 C collisions at 1.9A GeV with experimental collisions at 2.09 A GeV with experimental data, 

data (ref. 22). Dotted curve is proton exchange 
term with FSI, dash-dot curve is triton exchange 
with FSI, dashed-line curve is full calculation ne- 
glecting interference effects, and solid curve includes 
interference. 

The calculations in figures 20-23 show that a PWIA extraction of the internal momentum 
distribution is not possible in inclusive a-nucleus scattering. However, the agreement achieved 
suggests that a reasonable wave function has been used because several scattering mechanisms 
contribute and provide a strong constraint on model wave functions. The overlap function 
employed in our calculations was also found to give good agreement to experiments for pion- 
induced breakup at 5 GeV/c in reference 48. We also can conclude that a wave function with 
a minimum below 0.5 GeV/c is not in agreement with the a— 12 C data (refs. 53 and 58). 
The corrections to the impulse diagrams considered will not hide such a minimum and do 
not appear in the data. In figure 24, we compare the momentum distribution used here with 
those extracted in coincidence measurements with 426 MeV electrons (ref. 59) and 500 MeV 
protons (ref. 21). These experimental distributions will contain distortion effects peculiar to 
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Figure 22. Comparison of calculations of transverse 
momentum distribution for 3 H production in or- 12 C 
collisions at 0.385^4 GeV with experimental data. 


Figure 23. Comparison of calculations of transverse 
momentum distribution for 3 H production in a— 12 C 
collisions at 1.041 A GeV with experimental data. 



Figure 24. Momentum distribution for 3 H in 4 He 
ground state from ( e , ep) and (p, 2 p) reactions com- 
pared with distribution used in calculations. 


Total 


First order 

Secondrorder inelastic 



p, GeV/c 


Figure 25. Double-cluster scattering contributions to 
3 He production on 12 C at 385^4 MeV. 


the kinematics employed. Nevertheless, the comparison suggests a common shape out to about 
0 4 GeV/c; 

The contributions from the double-cluster scattering type of terms (eqs. (109) (112)) are 
illustrated in figures 25-26 for 3 He production on 12 C. The contribution from these rescattering 
c orre ctions is seen to be small in the longitudinal distribution; however they should increase, in 
importance at larger angles and for lower energies. 

The momentum distribution for a deuteron pair in ^He was fit to the Monte Carlo calculations 
of reference 10 as shown in figure 27. Calculations are shown in figures 28 and 29 with 
contributions from the various impulse terms shown. A significant contribution is seen for 
elastic fragmentation at forward angles. 


Total 




Figure 26. Double-cluster scattering contributions to 
3 He production on 12 C at 1041,4 MeV. 


Figure 27. Comparison of momentum distribution for 
d-d in 4 He Monte Carlo calculations of reference 10 
with three-term Yukawa fit. 


Jotal 

4 He — > dd elastic 

4 He — > dd inelastic 

4 He — » dnp 

4 Experimental with error 
(refs. 55 and 56) 



Total 



p, GeV/r 


Figure 28. Inclusive d production in a-\-p reaction at 
1.05,4 GeV. 


Figure 29. Inclusive d production in c*- 12 C collisions at 
2.09,4 GeV. 


Overall, comparing our results with experimental data for 4 He fragmentation is encouraging. 
We can identify several areas that will lead to greater predictive capability. The inclusive 
scattering model for cluster-target interactions should be extended to include charge exchange 
which will become important below several hundred MeV/amu. Also, these distributions should 
be extended to include pion production in order to make predictions at higher energies. Pion 
absorption on alpha clusters (^He, ^H, and 2 H) should then be studied. Final-state interaction 
effects in 2 H production must be included in order to make predictions of total-production 
cross sections. A first attempt will consider just FSI between two fragments while summing 
all contributions. Also, off-shell effects in higher partial waves and spin-orbit coupling effects 
on the FSI should be estimated. The two-body amplitudes employed should be improved for 
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accuracy at large momentum transfers, especially for 1 H targets, and should include nuclear 
medium effects that will become important at lower energies. Finally, an improved treatment of 
the many-body phase space for inelastic fragmentation would result from the use of equation (93) 
with the distribution of equation (75) rather than the approximation of the three-body phase 
space (eq. (94)) employed here. 


5. Parametric Data Base 

5.1. Interaction Cross Sections 

We next discuss parameterizations of interaction cross sections and energy spectra for alpha 
particles in common shielding materials. The work of Meyer (ref. 2) gives a complete summary 
of a- 1 !! cross sections based on measurements up to 1972. Parameterizations of H and He 
production on ! H below 300 MeV/amu were discussed in reference 60. An extensive list of 
earlier references of experiments is given in reference 2. More recent experiments are absorption 
cross sections between 18 and 48 MeV in reference 61. In reference 62, deuteron production 
at IAA GeV was measured with the inclusive deuteron production cross section reported at 
30.64 ± 0.62 mb. Also, postdating the compilation by Meyer is the result for A = 3 fragments in 
references 49 and 50. The most important shortcomings of the data base for He- H interactions 
are high-energy measurements above a few GeV and a complete absence of data for nucleon 
production cross sections. 

By using our theoretical estimates and the existing data, we parameterize the fragmentation 
cross sections for 3 He, 3 H, and 2 H production on 'Has 


asjje = 42.5 


h + e (T th -T)/6.8 
x ( 1 + 0.36 


1 ( 051 Y 

l J V 1 + 6.7e _r / 34 / 


T \ — (T- 780)/2300 

520 ) 


(127) 




= 15.5 


[ l + e (T th -T)/7 


-*](*" 


0.45 y 
+ 7e _r / 55 / 


x I 1 + I' 8 !/ 5^0 ) e- (T ~ 750)/4500 


(128) 


and 


Vd = 17 


Li + e (%-r)/i2 


- 1 


, . 0.21[(T/145) - 1] \ __r/; 3000 

A T" fiACL^/a I e 


1 + e (145-T)/ 6 


(129) 


where 7)/, is the threshold energy for the breakup reaction listed in table 5 (ref. 63) and T is the 
kinetic energy in units of A MeV. The low-energy behavior of equation (127) resembles that of 
reference 60. The pickup cross section is parameterized as 


^PICKUP 


= 48e“ 


(T-T^) 1 ' 7 / 1350 


(130) 


and contributes to both the inclusive 3 He and 2 H production cross sections. At low energies 
the resonance 5 Li occurs which is not considered here. The parameterizations are compared 
with experimental data in figures 30 32. AH cross sections are set constant above 3A GeV. The 
energy variations near thresholds and the pion production region are accurately reproduced. 



Table 5. Thresholds and Q - Values for p+ 4 He 


Reaction 

Q, MeV 
(a) 

Threshold, MeV 

4 He(p, df He 

-18.354 

22.94 

4 He(p, 2p) 3 He 

-19.815 

24.77 

4 He(p,pn) 3 He 

-20.578 

25.72 

4 He(p,pd) 2 H 

-23.848 

29.81 

4 He(p,ppn) 2 H 

-26.072 

32.59 

4 He(p,pprm) I H 

-28.297 

35.37 


a Data taken from references 60 and 63 . 


Total 

Stripping 




Figure 30. Comparison of parametric model for 3 He Figure 31. Comparison of parametric model for 3 H pro- 
production in a—p collisions with experimental duction in a— p collisions with experimental data, 

data. 

Total 

Stripping cross section 



Figure 32. Comparison of parametric model for 2 H production in a~p collisions with experimental data. 
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For the J H target the absorption (ABS) cross section below 80 MeV is assumed as 

i^ABS = ^3 He + ^3 H + + ^PICKUP ( 131 ) 

Above 80 MeV (and below 80 MeV for A T > 1), the energy-dependent parameterization of 
Townsend and Wilson (ref. 64) is used which is given by 


o-ABS = 107r 0 ( E ) 


■^Re 1-26 



(132) 


where 

0(E) = 1 + | 

6(E) = 0.2 + -h + y 0.292e -T / 792 cos(0.229T° 453 ) 


(133) 


with a normalization correction of 0.95 used for and the argument of the cosine function given 
in radians. In equation (132) the nuclear matter radii are used. The absorption cross section for 
is shown in figure 33 which presents an excellent reproduction of the experimental data. 

The proton and neutron production are expected to rise dramatically above pion production 
thresholds since two-body collisions will be predominantly inelastic leading to pion absorption 
in the A = 3 or A = 2 clusters. Pion production has not been treated in our theoretical 
considerations. However, the expected rise in proton and neutron cross sections occurs if we 
simply balance the absorption cross section with channels that do not lead to proton and neutron 
production, respectively, as shown in figure 34. 


Calculations 



Figure 33. Comparison of parametric model for a+ 1 H 
absorption cross section with experimental data. 



Figure 34. Parametric model predictions for proton 
production and neutron production in a+*H 
collisions. 


The experimental data base for composite targets is extremely small. A comparison of 
equation (132) with data for a+ 12 C absorption cross sections is shown in figure 35. Agreement 
is excellent and previous analyses (ref. 63) suggest similar agreement for other targets. The 
stripping reactions become more complicated for At >1 because several of these channels exist, 
and stripping or pickup to excited states of the target contributes to the complication. We 
follow Serber (ref. 25) and assume a surface reaction for nucleon stripping on He and then 
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Figure 35. Comparison of parametric model for a+ 12 C 
absorption cross sections with experimental data. 



Figure 36. Comparison of fits for 3 H and 3 He produc- 
tion cross sections with experimental data. 


1 / Q 

equation (130) is scaled by Ayi . Since 4 He is its own mirror nuclei, we will ignore coulomb 
effects and assume that equation (130) is used for both 3 He and 3 H production in stripping 
reactions. Note that 3 H is not produced in stripping on 1 H. A slight overestimate may occur 
because a small contribution from 3 He exchange is expected in the reaction — > 3 He+d. 

We also ignore any d production in the stripping process for Aj > 1. In figure 36 the results for 
3 H and 3 He production on 12 C are shown. The fragmentation cross section is scaled as A^ 31 for 
these fragments. The measurements are from Webber (ref. 50), and the data at 3.6 A GeV are 
from reference 65. 

In table 6 we compare parametric fits to secondary yields for charge fragments at 3. 6 A GeV 
for several targets as measured in reference 65. The experiment of reference 65 measured 
only peripheral events with detection angles less than 5°. We expect the measurements of 
1 H secondaries to be underestimates. The multiplicity for nucleon production from 4 He at high 
energies is between 1 and 1.2 as compared with a value of 2 assumed in existing cosmic ray 
codes. For 2 H we have used a scaling of A^ 4 from our parameterization in equation (129). 


Table 6. Comparison of Experimental Fragmentation Cross Sections for 4 He 
Projectiles at 3.6.4 GeV With Model Fits 

[Calculations are given in parentheses] 


Fragment 

Values of cr/r, mb, for targets 3 of — | 

Li 

c 

A1 

Cu 

l H 

2 H 

3 h 

3 He 

166 ± 13 (536.6) 
84 ± 15 (68.2) 

47 ±5 (52.7) 

48 ±5 (48.1) 

227 ±20 (592.0) 
91 ± 27 (91.2) 

58 ± 9 (65.4) 

49 ± 8 (59.6) 

319 ±34 (823.9) 
113 ±38 (128.2) 
73 ± 20 (84.1) 

70 ± 15 (76.7) 

417 ±45 (1294.9) 
159 ±45 (184.2) 

95 ±14 (109.9) 

95 ± 20 (100.2) 


“Data were taken from reference 65 which measured only particles in peripheral events (# < 5°). 
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5.2. Secondary Energy Spectra 

The energy spectra of secondary particles is expected to be governed by the internal 
momentum distribution of the projectile nucleons with a weak dependence on target mass 
number from dynamic effects. At low energies, kinematic restrictions limit the energy losses 
possible beyond any restriction provided by the internal momentum distribution. The energy 
spectra of secondaries from fragmentation is parameterized as 

= exp f-(T - T 0 + 6 F ) 2 /2Wl 1 (134) 

dEp V 7T Wp L J 

where the width W F , downshift S F , and beam energy T 0 are in units of A MeV. Values derived 
from integrating the inclusive momentum distribution of equation (112) over all angles are given 
in table 7 for several energies and targets. Comparisons of the fit provided by equation (134) 
with calculations are shown in figures 37-40 and are quite accurate. Note that the internal 
motion of the projectile constituents leads to fragments produced with velocities higher than the 

beam velocity. 

Table 7. Spectrum Parameters for 4 He Fragments 


(a) A = 3 fragments (b) A — 1 fragment 



The energy spectrum of the elastically scattered alpha particles is parameterized using the 
Bom term of the optical model expansion that is normalized to the coherent model results 
(ref. 66). A similar approach is followed to parameterize the quasi-elastic energy spectrum. 
Assuming a Gaussian density matrix for the target gives (from eq. (76)) 

= 2mjy2<TiN exp [-2m^z (u> - e^)] 0 (u> - e jgJ (135) 
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Figure 37. Comparison of calculation of energy spec- 
trum for 4 He+ 12 C — ► 3 H+X at 520 A MeV with 
Gaussian fit of equation (134). 



Figure 38. Comparison of calculation of energy spec- 
trum for 4 He+ 16 0 — ► 3 H+X at 1000 A MeV with 
Gaussian fit of equation (134). 



E, MeV/amu 

Figure 39. Comparison of calculation of energy spec- 
trum for 4 He+ 27 Al -► 3 H4-V at 1000A MeV with 
Gaussian fit of equation (134). 



E , MeV/amu 

Figure 40. Comparison of calculation of energy spec- 
trum for 4 He+ 1 H — > 2 H+V at 500 A MeV with 
Gaussian fit of equation (134). 


with R a and R t denoting the matter radii of the alpha particle and target, respectively, and 
B denoting the slope parameter. Equation (135) is expected to underestimate the spectrum 
at large values of uj because of multiple scattering and perhaps pion production. In figures 41 
and 42, illustrations of the fit are shown. 
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(0, MeV/amu (0, MeV/amu 

Figure 41. Parametric fits to a+ 27 Al -* a+X reaction Figure 42. Parametric fits to a+ 16 0 -» a+X reaction 
at 100, 300, 1000 A MeV. at 100, 300, and 1000A MeV. 

6. Concluding Remarks * 

An overview has been presented of calculations of interaction cross sections for high-energy 
alpha particles colliding on nuclear targets. Models of elastic, quasi-elastic, and fragmentation I 

channels were described using multiple scattering theory in the impulse approximation. A § 

discussion of future theoretical emphasis was given with pion production at high energies : 

and nuclear medium effects at low energies identified as the principal areas of future work. . 

Extensive comparisons were made with existing experimental data for high-energy alpha particles | 

interacting with nuclear targets, and good agreement was found. Parametric energy-dependent 
interaction cross sections and energy spectra are presented and discussed. 

NASA Langley Research Center 
Hampton, VA 23681-0001 
December 11, 1992 
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